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Abstract. A very important property of a statistical dis- 
tribution is to know whether it obeys Gaussian statistics 
or not. On the one hand, it is of paramount importance in 
the context of CMB anisotropy studies, since deviations 
from a Gaussian distribution could indicate the presence of 
uncorrected measurement systematics, of remaining fluc- 
tuations contributed by foregrounds, of deviations from 
the simplest models of inflation or of topological defects. 
On the other hand, looking for a non-Gaussian signal is 
a very ill-defined task and performances of various assess- 
ment methods may differ widely when applied in different 
contexts. 

In previous studies, we introduced a sensitive wavelet- 
based method which we apply here to the COBE/DMR 
data set, already extensively studied, using different ap- 
proaches. This provides an objective way to compare 
methods. It turns out that our multi-scale wavelet de- 
composition is a sensitive method. Yet we show that the 
results are rather sensitive to the choice of both the de- 
composition scheme and the wavelet basis. We find that 
the detection of the non-Gaussian signature is "marginal" 
with a probabihty of at most 99%. 
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1. Introduction 

The COBE-DMR data (Smoot et al. 1992) give the dis- 
tribution of temperature anisotropics of the Cosmic Mi- 
crowave Background (CMB) on most of the sky with a 
resolution of about 10°. This data set was the first to put 
constraints on the spectral index and the amplitude of the 
initial power spectrum (Fixsen et al. 1996). Further anal- 
yses of the CMB anisotropics should also give hints on the 
dominant mechanism which generates the initial density 
perturbation. As a matter of fact, the simplest inflation 
models generally predict a normal (Gaussian) distribu- 
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tion of the perturbations whereas the topological defects 
generate a non-Gaussian distribution. 

Many studies aimed at characterising the statistical 
properties of the COBE/DMR data (Amendola 1996; 
Kogut et al 1996; Ferreira et al. 1998; Heavens 1998; 
Pando et al. 1998b; Schmalzing & Gorski 1998; Brom- 
ley & Tegmark 1999; Magueijo 2000; Mukherjee et al. 
2000) . These investigations have relied on various statisti- 
cal indicators both in the real space (e.g. high order mo- 
ments of the temperature distribution, Minkowski func- 
tionals, genus analysis, ...) and in the dual, Fourier or 
wavelet transformed, space (e.g. N-point correlation func- 
tions, cumulants and moments of coefficients, ...). The 
COBE/DMR temperature fluctuations have long been 
considered Gaussian distributed. However more recently, 
their statistical nature has been very debated, as some of 
the groups mentioned above found a non-Gaussian signa- 
ture in the COBE/DMR four year maps. 

Regardless of the implications of such a signature, we 
use the COBE/DMR data as the only available test bench 
for applying, on real maps, a method proposed in Aghanim 
& Forni (1999) and Forni & Aghanim (1999), geared at 
detecting non-Gaussian features. This method is based on 
the statistical analysis of the wavelet coefficients obtained 
after a multi-scale decomposition of the signal. In fact, 
the non-Gaussian signatures are enhanced in the wavelet 
space, which makes an analysis in this space more sensi- 
tive to the statistical signatures. In this context, the major 
differences between the methods based on wavelets come 
from the choices of the decomposition scheme and the 
wavelet basis. We briefly describe these choices and the 
method in section 2. We present our results and discuss 
them in respectively section 3. and section 4. Finally, we 
give our main conclusions in section 4. 

2. Data and analysis method 

The COBE/DMR four year data was "cleaned" from the 
contribution of the Galactic foreground emissions using 
two methods. One is based on a combination technique 
(linear combination of the DMR maps) and results on 
what we refer to as the DCMB maps. The second tech- 
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nique is based on the subtraction of the Galaxy contribu- 
tion and it gives the DSMB maps. 

Wavelet analyses were shown to be very useful tools to 
detect non-Gaussian signatures on simulated CMB maps 
(Aghanim & Forni 1999; Hobson et al. 1999). They were 
apphed on the COBE/DMR maps by Pando et al. (1998b) 
and more recently in an exhaustive study by Mukherjee 
et al. (2000). In this paper, we re-analyse both sets of 
COBE/DMR data (DCMB and DSMB) to study their 
statistical nature (Gaussian or non-Gaussian). In order for 
our conclusions to be the most independent possible of the 
cleaning method, we have analysed the less contaminated 
regions of the COBE/DMR sky namely: the northern and 
the southern poles. These regions arc associated respec- 
tively with the faces #1 and #6 of the so-called COBE 
cube on which the data were projected. 

Based on the method proposed in Forni & Aghanim 
(1999), we decompose the signal on the wavelet basis and 
compute the kurtosis which is the fourth order moment 
of the wavelet coefficient distribution. We thus derive the 
excess of kurtosis that is the deviation from the predicted 
value in a Gaussian distribution. The values obtained 
are compared with the excesses of kurtosis calculated us- 
ing 2000 Gaussian realisations of the COBE/DMR maps. 
These Gaussian realisations have the same power spectra 
as the COBE/DMR maps with truly random phases, ev- 
erything else (pixel sizes, coverage) being kept the same. 
Any statistical difference between the excess of kurtosis 
measured on COBE/DMR and those measured on the as- 
sociated Gaussian realisations should thus be associated 
with the non-Gaussian character of the data. 

The wavelet analysis consists in hierarchically decom- 
posing an input signal into a series of successively lower 
resolution "reference" signals and their associated "detail" 
signals (Mallat 1989). At each decomposition level, L, the 
reference signal has a resolution reduced by a factor of 
2^ with respect to the input signal. Due to the orthog- 
onality properties of the basis, a function can be com- 
pletely characterised by the wavelet basis and the wavelet 
coefficients of the decomposition. We test the detection of 
non-gaussianity against the wavelet basis in the "dyadic" 
decomposition scheme. We will compare our results with a 
study based on an alternative decomposition scheme: the 
"pyramidal" decomposition (Mukherjee et al. 2000). The 
dyadic decomposition is a transform in which, at each level 
L, only the reference signal (low-pass part of the signal) 
is decomposed. At each decomposition level, the analy- 
sis is applied in both directions and the total number of 
sub-bands after L levels of decomposition is 3L-I-1. The 
pyramidal decomposition is similar to the previous decom- 
position in the sense that only the reference sub-band is 
decomposed at each level. However, it is performed here 
separately in the two directions of the image. The total 
number of sub-bands after L levels of decomposition is 
then (L-hl)2. 
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Table 1. For the north galactic pole and the Haar wavelet 
basis, the first column represents the mean excess of kurto- 
sis computed with the wavelet coefficients over 2000 Gaus- 
sian reahsations of the COBE/DMR #1 map (DCMB). 
(T_ and (7+ are the standard deviations with respect to 
the mean of the kurtosis distribution. The fourth col- 
umn represents the excess of kurtosis computed on the 
COBE/DMR map and the last column gives the proba- 
bility, computed with the corresponding distribution, for 
the excess of kurtosis kQOBE ^'^ indicate a non-Gaussian 
signal. 



The choice of the wavelet basis is critical to test for 
non-gaussianity. We expect that the non-Gaussian fea- 
tures are due to point sources or sharp gradients in the 
signal. The detection of non-gaussianity is based on the 
search of these gradients. We therefore choose to use anti- 
symmetric wavelets which are proportional to, or almost 
proportional to, a first derivative operator and thus par- 
ticularly sensitive to gradients. In practice, we choose the 
6/10 tap filter (filter #3) of Villasenor et al. (1995) mainly 
used in data compression. It has, in fact, a very good im- 
pulse response and a low shift variance which better pre- 
serve the amplitude and the location of the details. We 
have also used the commonly used Haar wavelet which is 
also an anti-symmetric function. 

In the dyadic wavelet decomposition and using the 
anti-symmetric wavelet basis, we can discriminate between 
the coefficients associated with vertical, horizontal and di- 
agonal gradients. They are analogous to the partial deriva- 
tives, respectively d/dx, d/dy and d^ /dxdy of the signal. 

The maps we analyse have 32 x 32 pixels, with a pixel 
of 2.8° aside. Due to this reduced number, we restrict to a 
three level decomposition. The three decomposition scales 
(hereafter I, II and HI) correspond respectively to an an- 
gular scale 9 — 5.6° (that is a multipole order i = 35.7), 
e = 11.2° {I = 17.8), and 9 = 22.4° {I == 8.9). 



3. Results 
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Fig. 1. For the north galactic pole of the COBE/DMR #1 map 
(DCMB) and the Haar wavelet basis, we show the distribution 
of the excesses of kurtosis computed with the wavelet coeffi- 
cients (d/dy) at the third decompostion scale for 2000 Gaus- 
sian realisations. The horizontal dashed line represents the ex- 
cess of kurtosis measured on the COBE map with a probability 
of 98.7% (cf. Ime 8, Table^. 



3.1. North galactic pole 

We test for the detection of a non-Gaussian signature in 
the COBE/DMR, DCMB map, of the north galactic pole 
against 2000 Gaussian realisations of this map. First using 
the Haar wavelet basis in a three level decomposition, we 
compute the excesses of kurtosis of the coefficients asso- 
ciated with the details for the COBE map {k^Q^^, see 
Table |l|), and the Gaussian reahsations. We give, in this 
case, the mean excess of kurtosis k and the standard devi- 
ations (T_ and (7+ with respect to this mean for the quan- 
tities that are respectively smaller and larger than k. Fi- 
nally, we compare the obtained results and quantify the 
probability P that the measured excess of kurtosis on the 
COBE map [k^Qg^) indicates a non-Gaussian process. 
This probability is computed by integrating the excess of 
kurtosis distribution for the Gaussian realisations, one ex- 
ample of which is shown in Figure yj. It ranges between 
about 50 and 99% depending on the decomposition scale 
and on the type of details we study. We note that the non- 
Gaussian signature is detected at all scales with at least 
a probability larger than 68%. However, the most signifi- 
cant detections at about the 98 and 99% confidence levels 
are obtained respectively at the second and third scales 
where the beam and noise effects are expected to be less 
dominant. In the case of the Villasenor et al. wavelet ba- 
sis, the probability of detecting a non-Gaussian signature 
in the DCMB COBE/DMR #1 map lies between about 
67 and 98% (see Tabic pi). Again this detection is made at 
all scales. 

We perform the same analysis using the two wavelet 
basis on the COBE map obtained with the second sub- 
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Table 2. Same as Tabic 0, for the the Villasenor et al. 
wavelet basis. 



traction technique (DSMB maps). In Table |^, we give the 
probability for the excess of kurtosis, measured on the 
COBE/DMR map, to be associated with the detection of 
a non-Gaussian signal. We note again that the probabil- 
ities indicate a non-Gaussian process at all scales. The 
Haar basis seems however to give smaller (less discrimi- 
nating) values than the Villasenor et al. basis. 

3.2. South galactic pole 

Using the power spectrum of the COBE/DMR fj^Q map 
(DCMB) of the south galactic pole, we again simulate 
2000 Gaussian realisations of this map. We perform the 
wavelet decomposition using the Haar and the Villasenor 
et al. bases, and analyse the statistical properties of each 
type of details. As for the north pole, we compute the 
excess of kurtosis in the COBE/DMR map (k^Q^^) and 
the probability, P, that it belongs to the distribution of 
the excesses of kurtosis of the Gaussian realisations. The 
latter is characterised by its mean k and standard devia- 
tions CT- and (7+ of respectively smaller and larger values 
with respect to the mean. The results are summarized in 
Tables | and |. 

At first sight, these results suggest that the southern 
hemisphere distribution is more Gaussian than its north- 
ern hemisphere counterpart. As a matter of fact, the prob- 
ability that the excess of kurtosis ^cobb measured on the 
COBE DCMB map indicates a non-Gaussian process lies 
between 53 and 92% for the Haar wavelet basis and be- 
tween 47 and 87% for the Villasenor et al. wavelet. At all 
scales, we measure a probability larger than 68% for at 
least one type of details. However, these values are always 
smaller than those obtained for the northern hemisphere. 

When we analyse the DSMB map of the southern 
hemisphere (Table ra), we note the same behaviour: the 
non-Gaussian signature is weaker in the south than in the 
north. The probabilities lie between 44 and 87% for the 
Haar wavelet basis, and between 45 and 98% for the Vil- 
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Table 5. Probability for the excess of kurtosis of the COBE/DMR map (DSMB) to be attributed to a non-Gaussian 
signature. The numbers are given, for the three types of details, as a function of the decomposition scale and the 
wavelet basis. The first part of the table represents the north pole and the second part is for the south pole. 



Scale I 


k 


0"_ 


0-+ 


f^COBE 


P 


d/dx 

d/dy 

d"^ /dxdy 


0.006 
-0.01 
0.02 


0.26 
0.26 
0.27 


0.35 
0.36 
0.35 


-0.28 
-0.21 
-0.28 


85.7% 
75.8% 
85.7% 


Scale II 


1 


d/dx 

d/dy 

d"^ /dxdy 


0.03 
0.003 
-0.03 


0.50 
0.49 
0.46 


0.82 
0.80 
0.75 


-0.56 
0.59 
-0.62 


87.6% 
87.7% 
91.5% 


Scale III 


1 


d/dx 

d/dy 

d^ /dxdy 


-0.03 
0.007 
0.005 


1.01 
0.99 
1.03 


1.89 
1.81 
1.93 


-0.68 
-1.27 
0.40 


67.5% 
90.3% 
52.9% 



Table 3. For the south galactic pole and the Haar wavelet 
basis, the first column represents the mean excess of kurto- 
sis computed with the wavelet coefficients for 2000 Gaus- 
sian realisations of the COBE/DMR #6 map (DCMB). 
(T_ and (T+ are the standard deviations with respect to 
the mean of the kurtosis distribution. The fourth col- 
umn represents the excess of kurtosis computed on the 
COBE/DMR map and the last column gives the proba- 
bility, computed with the corresponding distribution, for 
the excess of kurtosis k^oBE ^^ be associated with a non- 
Gaussian signal. 
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Table 4. Same as Table |^ for the Villasenor et al. wavelet 
basis. 



lasenor et al. basis. The highest probabilities, in this case, 
are obtained at the second decomposition scale. 

3.3. Scale-scale correlation spectra 

Another method to detect a non-Gaussian signature 
consists in using the scale-scale correlations as pro- 
posed by Pando et al. (1998a), Pando et al. (1998b) and 



Mukherjee et al.(20GG). The correlations between struc- 
tures at adjacent decomposition scales are quantified. 
They can be computed using the expression given in 



Mukherjee et al.(2000) 
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In this expression, the hj^j^^i^^i^ represent the coefficients 
in the wavelet decomposition. Cj^j^ gives the correla- 
tion factor between the coefficients in the domain {ji,J2) 
and its adjacent scale {ki,k2), where fci = ji + 1 and 
^2 — ^2 + 1, for a Gaussian signal Cj-^j^ = 1. In or- 
der to assess the gain brought by this further step of 
analysis, we proceed as previously. We first compute the 
scale-scale correlation coefficients for our 2000 Gaussian 
realisations of the COBE/DMR maps after decomposing 
them on three scales with the Haar and Villasenor et al. 
bases. We then compute the scale-scale correlations for 
the "real" COBE maps and finally compute the probabil- 
ity that the measured numbers indicate a non-Gaussian 
signature by using the distribution of correlations for the 
Gaussian maps. 

The results are summariz ed in Table p]. They bas ically 
agree with those obtained by Mukherjee et al.(2000) when 
the comparison can be made, that is for the diagonal de- 
tails. The probabilities for the correlations between scales 
I and II in the north galactic pole lie between 64 and 76%, 
with the Haar basis, and between 62 and 99% with the Vil- 
lasenor et al. basis. When we correlate scales II and HI, 
we find a probability of 69 to 96%, with Haar; whereas 
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Table 6. Scale-scale correlation coefficients for the three types of details as a function of the wavelet basis. For the 
south and north pole, we give the probability for the correlation coefficient to be attributed to a non-Gaussian signal. 
The correlations are performed between the first and second scales, and between the second and third scales. The 
upper part of the table is given for the DSMB COBE/DMR maps whereas the lower part is for the DCMB maps. 



it is of 61 to 88% with the Villasenor et al. basis. In the 
south galactic pole, the probability for the correlation co- 
efficients between scales I and II to be associated with a 
non-Gaussian signal ranges between 60 and 94% with the 
Haar basis, and between 63 and 98% with the Villasenor 
et al. basis. Correlating scale II and scale III, the proba- 
bilities are between 44 and 87% with Haar and between 50 
and 97% with Villasenor et al. The results obtained with 
the correlation coefficients agree with the direct analysis of 
the wavelet coefficients. However, they seem to give larger 
non-Gaussian signatures for the south pole than the direct 
analysis. 

4. Discussion and conclusions 

At each decomposition scale, we detected a non-Gaussian 
signature using three different criteria related to the hor- 
izontal, vertical and diagonal details of the image. We 
tested and compared the sensitivity to non-Gaussian fea- 
tures of the south and north pole maps obtained after sub- 
tracting the Galaxy with two techniques. We thus anal- 
ysed this four map set with two wavelet bases used in a 
three level dyadic decomposition. All the results are sum- 
marized in Tables 0, |, |, |, |. 

This rather exhaustive study leads to two main con- 
clusions. One is that the analysis of the COBE maps is 
rather sensitive to the choice of the wavelet basis when 
the same decomposition scheme is used. The second con- 
clusion is that the non-Gaussian signature is detected 
with a probability of 99%. At the first and second de- 
composition scales, the high probalities are very likely 
associated with some instrumental effects, such as pixel 
to pixel correlations or Poisson noise, or systematic ef- 
fects (Banday et al. 2000). At the third decomposition 
scale (which is less affected by the instrumental effects), 
we detect a non- Gaussian signature with a probability < 
95% in the southern hemisphere, whereas in the north- 
ern hemisphere, it is detected with a probability close to 



99%. We have compared our results to those obtained by 



Mukherjee et al.(2000). Namely, we compare the results 
they obtain for what they note k — 32, 16 and 8, which 
correspond exactly to what we refer to as the diagonal 
details for respectively the first, second and third decom- 
position scales. We used the dyadic decomposition advo- 
catcd by [Aghanim et al.(1998| |Forni fc Aghanim(1999) 



Aghanim & Forni(1999) r ather than the pyramidal 
decomposition used in [Hobson et al.(1999) and in 
Mukherjee et al.(2000) . This choice is motivated by the 
fact that this particular decomposition scheme gives the 
correlations between the two directions of the map at each 
scale. It is not the case with the pyramidal decomposition 
which treats both directions as if they were independent. 
Another advantage of the dyadic decomposition is that, 
for each decomposition scale, the analysis bears upon the 
maximum number of coefficients possible. This is crucial 
for statistics in general and for the small COBE maps (#1 
and 6) in particular. 

The comparison between our results, obtained with 
the Villasenor et al. wavelet basis, and those of 



Mukherjee et al.(2000) shows that our non-gaussianity de- 



tection technique gives every time higher confidence lev- 
els except for two cases which are the first decomposi- 
tion scale of the southern hemisphere (DCMB and DSMB 
maps). When we use the Haar wavelet basis, our confi- 
dence levels are higher than, or similar to, those quoted in 
Mukherjee et al.(2000) except for two cases that are the 
third decomposition scale of the northern and southern 
hemisphere of the DSMB map. Of course in addition to 
the diagonal details, we have two other criteria to detect 
the non-Gaussian signature (vertical and horizontal de- 
tails). In particular, for the cases where the probability 
is smaller than in Mukherjee et al.(2000) we find a non- 
Gaussian signature based on the vertical or/and the hor- 
izontal details. 

In our study, we find that the Villasenor et al. wavelet 
basis combined with a dyadic decomposition gives higher 
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probabilities than the pyramidal decomposition. This ba- 
sis is also more sensitive than the Haar basis. The better 
performances in detecting gradients come from the larger 
size of its kernel (Villasenor et al. 1995). However, this 
could become a weakness when a small number of pixels 
is available (small maps or high decomposition levels). In 
this case, the Haar wavelet basis could be more adapted 
due to the smaller size of its kernel. 

Our choice of wavelet basis and decomposition scheme 
confirms that a non-Gaussian signature can be detected ei- 
ther in the south or the north galactic pole. It seems how- 
ever that the statistical signature is less significant in the 
southern pole. In this case, the largest probability, > 95%, 
is obtained only at the second decomposition scale for the 
DSMB COBE/DMR #6 map. The non-Gaussian signa- 
ture is more significant for the north galactic pole as the 
probability is often > 95% and can even reach 99%. The 
scale-scale correlations represent a step further in the sta- 
tistical analysis of the signal. However in our case, it does 
not make a significant improvement in the results. The 
non-Gaussian signature is indeed detected with a proba- 
bility of 99% at most, such as in the direct analysis. 

The large size of the beam and pixels (respectively 7 
and 2. 8°) results in a rather small number of wavelet co- 



efScien lts. This, together with the low signal to noise ratio 
of the COBE/DMR data, represent very limiting factors 
which prevent us from drawing any hasty conclusion from 
the rather marginal detections of the non-Gaussian signa- 
tures. This kind of statistical analyses will soon demon- 
strate all their performances when applied to the new gen- 
eration of CMB data sets of higher sensitivity and resolu- 
tion (Boomerang (De Bernardis et al. 2000) and Maxima 
(Hanany et al. 2000)). We thus expect that these studies 
will then put constraints on the primordial non-Gaussian 
contributio ns such as the topolog ical defects suggested for 
example in Bouchet et al.(2000) . 
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